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Abstract 

Cellular and viral S-adenosylmethionine-dependent methyltransferases are involved in many regulated processes such as 
metabolism, detoxification, signal transduction, chromatin remodeling, nucleic acid processing, and mRNA capping. The 
Severe Acute Respiratory Syndrome coronavirus nsp16 protein is a S-adenosylmethionine-dependent (nucleoside-2'-0)- 
methyljtransferase only active in the presence of its activating partner nspJO. We report the nsp10/nsp16 complex structure 
at 2.0 A resolution, which shows nspIO bound to nsp16 through a ~930 A 2 surface area in nspIO. Functional assays identify 
key residues involved in nsp10/nsp16 association, and in RNA binding or catalysis, the latter likely through a SN2-like 
mechanism. We present two other crystal structures, the inhibitor Sinefungin bound in the S-adenosylmethionine binding 
pocket and the tighter complex nsp10(Y96F)/nsp16, providing the first structural insight into the regulation of RNA capping 
enzymes in (+)RNA viruses. 
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Introduction 

Most eukaryotic cellular and viral mRNAs are modified by the 
addition of a polyadenine tail at the 3'- terminal and a cap 
structure at the 5'-terminal. The RNA cap protects mRNA from 
degradation by 5' exoribonucleases, ensures efficient mRNA 
translation, and prevents recognition of viral RNA via innate 
immunity mechanisms [1,2,3,4], The RNA cap is made of an NT- 
methylated guanine nucleotide connected through a 5'-5' 
triphosphate bridge to the first transcribed nucleotide, generally 
an adenine. Through 2'-0 methylation of the latter, this cap-0 
structure ( 7Me GpppA...) may be converted into a cap-1 structure 
(7MeGpppA 2 , n Me ...). I n the eukaryotic cell, the cap is added co- 
transcriptionally in the nucleus by three sequential enzymatic 
reactions[1,5]: (i) an RNA triphosphatase (RTPase) removes the 5' 
“/-phosphate group of the nascent mRNA; (ii) a guanylyltransferase 
(GTase), dubbed capping enzyme, catalyses the attachment of 
GMP to the 5'-diphosphate mRNA; and (iii) an S-adenosylme¬ 
thionine (SAM)-dependent (N7-guanine)-methyltransferase 
(N7MTase) methylates the cap onto the N7-guanine, releasing 
S-adenosylhomocysteine (SAH). In general, a SAM-dependent 


(nucleoside-2'-0-)-methyltransferase (2'-0-MTase) further inter¬ 
venes, in higher eukaryotes, to yield a cap-1 structure. 

The viral RNA capping machinery is structurally and 
mechanistically diverse, and RNA viruses often deviate from the 
paradigmic eukaryotic mRNA capping scheme. For example, 
alphaviruses methylate GTP onto the N7-guanine before the 
presumed attachment of 7M ‘GMP to the nascent viral 5'- 
diphosphate mRNA [6]. In the case of single-stranded negative- 
sense (-)RNA viruses, such as the vesicular stomatitis virus, the L 
polymerase attaches GDP rather than GMP to a nascent viral 5'- 
monophosphate mRNA, covalently linked to the viral capping 
enzyme [7]. Other viruses, such as influenza virus capture a short 
capped RNA oligonucleotide from host cell mRNAs and use it as 
an RNA synthesis primer. This process is known as « cap 
snatching »[8]. 

In 2003, a novel coronavirus named Severe Acute Respiratory 
Syndrome coronavirus (SARS-CoV[9]) was responsible for the first 
viral pandemic of the new millennium with ~8000 cases globally 
and a 10 % case-fatality rate. Coronaviruses encode an unusually 
large membrane-associated RNA replication/transcription machin¬ 
ery comprising at least sixteen proteins (nspl-to-16)[10]. For 
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Author Summary 

A novel coronavirus emerged in 2003 and was identified as 
the etiological agent of the deadly disease called Severe 
Acute Respiratory Syndrome. This coronavirus replicates 
and transcribes its giant genome using sixteen non- 
structural proteins (nsp 1-16). Viral RNAs are capped to 
ensure stability, efficient translation, and evading the 
innate immunity system of the host cell. The nsp16 protein 
is a RNA cap modifying enzyme only active in the presence 
of its activating partner nsplO. We have crystallized the 
nsp10/16 complex and report its crystal structure at atomic 
resolution. NsplO binds to nsp16 through a ~930 A 2 
activation surface area in nsplO, and the resulting complex 
exhibits RNA cap (nucleoside-2'-0)-methyltransferase ac¬ 
tivity. We have performed mutational and functional 
assays to identify key residues involved in catalysis and/ 
or in RNA binding, and in the association of nsplO to 
nspl6. We present two additional crystal structures, that of 
the known inhibitor Sinefungin bound in the SAM binding 
pocket, and that of a tighter complex made of the mutant 
nsplO(Y96F) bound to nspl6. Our study provides a basis 
for antiviral drug design as well as the first structural 
insight into the regulation of RNA capping enzymes in 
(+)RNA viruses. 


SARS-CoV, the RNA cap structure likely corresponds to a cap-1 
type [11,12,13]. As in many other (+)RNA viruses, the RTPase 
activity is presumably embedded in the RNA helicase nsp 13, 
whereas the GTase remains elusive. RNA cap 2'-0-MTase activity 
was first discovered in the feline coronavirus (FCoV) nsp 16 [14], 
Shortly after, SARS-CoV nsp 14 was shown to methylate RNA caps 
in their N7-guanine position[15]. Curiously, although closely 
homologous to that of FCoV, recombinant SARS-CoV nsp 16 
alone was devoid of enzymatic activity. It was demonstrat- 
ed[16,17,18,19] that nsplO interacts with nspl6, conferring 2'-0- 
MTase activity to nsp 16 on N7-methyl guanine RNA caps 
selectively[16]. The latter selectivity implies that RNA cap 
methylation obeys an ordered sequence of events during which 
nspl4-mediated N7-guanine methylation precedes nspl0/nspl6 
RNA2'-0 methylation. NsplO is a double zinc finger protein of 148 
residues whose crystal structure is known[20,21]. Together with 
nsp4, nsp5, nsp 12, nsp 14, and nsp 16, nsplO has been found to be 
essential in the assembly of a functional replication/transcription 
complex[22]. Drawing on these observations, nsplO has been 
proposed to play pleiotropic roles in viral RNA synthesis [2 3] and 
polyprotein processing through interaction with the main protease 
nsp5[24], 

SAM-dependent MTases belong to a large class of enzymes 
present in all life forms. These enzymes catalyze the transfer of the 
SAM methyl group to a wide spectrum of methyl acceptors, 
indicating that a common chemical reaction is used on a variable 
active-site environment able to activate the methyl acceptor atom. 
Although SAM-dependent MTases share little sequence identity, 
2'O-MTases exhibit a KDKE catalytic tetrad and a very 
conserved folding made of a seven-stranded fl-sheet surrounded 
by one to three helices on each side [25], always similar to the 
paradigmatic catechol-0-MTase[26]. The SAM binding site 
general location is conserved, suggesting that evolutionary pressure 
on the MTase fold has maintained the same SAM-binding region 
whilst accommodating the versatile chemistry of the methyltrans- 
fer reaction. 

Structural and functional studies of viral MTases involved in 
RNA capping is an expanding research area, since these enzymes 
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show unexpected diversity relative to their cellular counterparts, 
and thus constitute attractive antiviral targets. Crystal structures of 
viral RNA cap MTases exist for only three viral families, namely 
Poxviridae , Reoviridae, and Flaviviridae. The Vaccinia virus VP39 
ciystal structure was the first to be elucidated in 1996 [2 7]. The 
structure of this DNA virus RNA 2'-0-MTase revealed a 
conserved MTase fold similar to that of RrmJ (also named FtsJ), 
the canonical reference folding for RNA cap MTases[26], More 
recently, the crystal structure of a second Vaccinia virus N7- 
guanine RNA cap MTase domain (Dl) was determined in 
complex with its activator protein D12[28], The study revealed 
that D12 also bears an MTase fold, but has lost catalytic capability 
due to truncation of its SAM binding site. In turn, Reoviridae 
provided the first RNA cap MTase structures at 3.6 A resolution 
as forming part of the reovirus core [29]. Another RNA cap 
machinery was more recently described for the non-turreted 
orbivirus Bluetongue virus VP4 protein at 2.5 A resolution[30], 
which revealed a three-domain protein, with a “head” guanylyl- 
transferase domain, a central N7-guanine MTase, and a “bottom” 
2'-0-MTase domain. This architecture illustrates the sequence of 
three out of the four chemical reactions involved in RNA capping 
described above. 

Regarding (+)RNA viruses, MTase structural information at the 
atomic level is only available for a single genus. The flavivirus N- 
terminus domain (residues 1-265) of the NS5 RNA-dependent 
RNA polymerase harbors an RrmJ fold with an N-terminus 
extension able to accommodate RNA cap structures [31,32]. This 
enzyme carries both N7-guanine MTase and 2'-0-MTase 
activities on a single domain with one shared active site [3 3]. 
Homologous domains have been crystallized for a number of 
flaviviruses, revealing a conserved fold and activity [34], suggesting 
that MTases might represent interesting targets for drug design. 
No other (+)RNA virus RNA cap MTase ciystal structures have as 
yet been defined. 

In 2003, the identification of the 2'-0-MTase signature 
sequence in the SARS-CoV genome added nsp 16 to the list of 
putative targets for antiviral drugs [35]. Several compounds have 
been shown to inhibit viral MTases, such as the co-product of the 
MTase reaction SAH, Sinefungin, and aurintricarboxylic acid 
(ATA)[14,36,37,38,39]. In this paper, we report the crystal 
structure of the SARS-CoV 2’-0-MTase nsp 16 in complex with 
its activator, the zinc finger protein nsplO, at 2.0 A resolution, in 
conjunction with mutagenesis experiments, binding and activity 
assays. These results lay down the structural basis for the nsplO 
function as an activator of nspl6-mediated 2'-0-MTase. We 
identify residues playing key roles in the nsp!0/nspl6 interaction, 
as well as other residues involved in 2’-0-MTase catalysis and 
RNA binding. We also report the crystal structure of the nsplO/ 
nsp 16 complex bound to the inhibitor Sinefungin. Comparison 
with known cellular SAM binding sites points to the nsp 16 
nucleobase binding pocket as a possible target for the design of 
selective antiviral molecules. 

Results 

Crystallization and Structure Determination of an Active 
nsp16 2'-0-MTase 

We observed that purified nsp 16 was unstable in solution, 
impeding crystallogenesis. Yeast double-hybrid and co-immuno- 
precipitation experiments on purified SARS-CoV nsplO and 
nsp 16 have uncovered the reciprocal interaction of these two 
proteins[16,17,18], Indeed, SARS-CoV nspl6 exhibits 2'-0- 
MTase activity only when complemented with SARS-CoV 
nsplO, raising the interesting possibility that nsplO acted as a 
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scaffold for nspl6. Co-expression of nsplO and nspl6 using a bi- 
cistronic prokaryotic expression vector facilitated affinity chroma¬ 
tography purification and crystallization of the complex[16,40]. 
Crystals diffracted to ~1.9 A. The position of the nsplO protein 
was determined using molecular replacement with the SARS-CoV 
nsplO protein structure[20] as a search model. Strong peaks in 
both the residual and anomalous Fourier maps confirmed the 
presence of two zinc ions. Nspl6 was well defined by its electron 
density except for two flexible loops (residues 19-35 and 135-137) 
with high B factors and weak or missing electron density. These 
loops are solvent-exposed at each side of the putative RNA- 
binding groove (see below). Structure determination data and 
refinement statistics are reported in Table 1. 

Structure of the nsp10/16 Heterodimer 

The heterodimer can be conveniently viewed as nsp 16 sitting on 
top of a nsplO monomer (Fig. 1A). The nsplO overall structure in 
the complex remains essentially unchanged relative to published 
structures of nsplO alone, with its N-terminus comprising two a- 
helices, a central (3-sheet domain, and a C-terminus domain 
containing various loops and helices (see[20,21], Fig. IB). 
Comparison with existing crystal structures of nsplO using 
DaliLite[41] rendered nsplO atomic coordinates very similar to 
those of nsplO in our nspl0/nspl6 complex. The average RMSD 
is about 0.77 A in 118 residues (PDB codes 2FYG, 2G9T and 
2GA6[20,21]). This indicates that neither significant conforma¬ 
tional change nor surface modification occurs in nsplO when 
binding to nsp 16. The nsplO structural Zn 2+ ions are not directly 
involved in the nspl0/nspl6 interface (Fig. 1A). 

Nsp 16 adopts a canonical SAM-MT fold (Figs. IB, 2A and B), 
as defined initially for the catechol 0-MTase[25]. The seven- 
stranded (3-sheet MTase fold has been described as having a 
secondary structure topology defining two binding domains, one 
for SAM and the other for the methyl acceptor substrate (Fig. 2A). 
The nsp 16 topology matches those of dengue virus NS5 


N-terminal domain and of vaccinia virus VP39 MTases[27,31]. 
Nsp 16 lacks several elements of the canonical MTase fold, such as 
helices B and C (Fig. 2B). 

S-adenosylhomocysteine- and Sinefungin-Binding 

Electron density corresponding to one molecule of S-adeno- 
sylhomocysteine (SAH), the co-product of the methylation 
reaction, was identified in the putative SAM-binding site 
(Figs. 1A and 3A). Neither SAM nor SAFf was added to the 
purification or crystallization buffers, therefore it must have been 
captured from the medium by nsp 16 during bacterial growth. The 
SAFI molecule is found with its adenine in an anti conformation 
and the ribose pucker in a southern (2’-endo/3'-exo) conformation. 
All the residues involved in SAM/SAH binding are absolutely 
conserved in coronavirus npl6s (Fig. SI). Binding specificity for 
SAM/SAH is achieved by holding distal SAM/SAH carboxylic 
and amino groups through five hydrogen bonds (G81, N43, Y47, 
G71, and D130) (Fig. 3A and Fig. S2A). The ribose moiety is held 
by three hydrogen bonds involving Y132, G73, and D99. As in the 
case of other MTases[25], the SAH binding cleft is globally 
positively charged. However, an aspartic acid (D99) acts as the 
ribose-sensing residue with its side chain carboxyl making strong 
hydrogen bonds with both ribose hydroxyls (Fig. S2A). Binding of 
the adenine base involves few contacts. The nucleobase occupies a 
loose hydrophobic pocket engaging two hydrogen bonds of 
moderate strength with side chain and main chain atoms of 
conserved residues D114 and Cl 15, respectively. 

Soaking the crystals into a Sinefungin-containing buffer 
captured this MTase inhibitor in the SAH binding site almost 
perfectly superimposable on SAH (Fig. 3B and Fig. S2B). Binding 
involved the same residues and contacts as SAH. Inhibition of the 
MTase reaction by Sinefungin dierefore probably occurs compet¬ 
itively. The Sinefungin amino group quasi-isosteric to the donated 
SAM methyl group indicates a cavity where the 2 '-hydroxyl of the 
capped RNA is expected to bind. Lining this empty substrate 


Table 1 . Crystal, collection, structure determination data and refinement statistics. 


DATA 

wild-type (SAH) 

nsp10(Y96F)/nsp16 

wild-type (SFG) 

Instrument 

SOLEIL, Proxima 1 

ESRF ID14-1 

ESRF ID23-1 

Wavelength 

0.9792 

0.9334 

0.9792 

Space group 

C222t 

C222t 

C222, 

Cell dimensions a, b, c (A) 

68.07 184.62 128.83 

68.15 184.80 129.01 

68.42 185.04 129.46 

Resolution range (A) 

37.52-2.00 (2.11 - 2.00)* 

45.51-2.05 (2.17 - 2.05) 

45.58 -2.50 (2.64 - 2.50) 

Total number of reflections 

201703 (29338) 

187379 (27241) 

211213 (30837) 

Number of unique reflections 

54947 (7963) 

51033 (7362) 

28921 (4163) 

Completeness (%) 

99.7 (100) 

99.6 (100) 

100 (100) 

1/aV) 

7.1 (2.7) 

7.5 (3.2) 

12.5 (4.3) 

Rsym ** 

0.114 (0.430) 

0.112 (0.350) 

0.102 (0.400) 

Multiplicity 

3.7 (3.7) 

3.7 (3.7) 

7.3 (6.7) 

Refinement (20 cycles) 

p*** 

0.213 

0.200 

0.201 

Rfree 

0.227 

0.234 

0.215 

RMSD bond length (A) 

0.007 

0.007 

0.011 

RMSD bond angle (A) 

1.015 

0.968 

1.250 


*Values in parentheses give the high resolution shell values. 
**Rsym =E |I-<I>|/Z I. 

***R = Z||Fo|-|Fc||/E |Fo|. 

doi: 10.1371/journal.ppat.1002059.t001 
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Figure 1. Cartoon representation of the nsp10/nsp16 complex with the reaction product SAH and metal ions. A) The nsp16 protein 
(green) is bound to nspIO (yellow) through an interface which does not involve Zn ions (grey spheres) present in nsplO. One metal ion (blue sphere) 
is found in nsp16 on the opposite face from the active site to which a SAH molecule is found (red sticks). B) Ribbon representation of nsp16, rainbow 
colors from N- to C-terminus. Top: Each secondary structure element is labeled according to[25] (see also Fig. S3). The SAH molecule shown in sticks 
colored following atom type. Bottom: View of interface involving nsp16 (green surface) and nsplO (yellow ribbons) showing the nspIO secondary 
structure elements involved in the interface. 
doi:10.1371/journal.ppat.1002059.g001 


cavity are the residues proposed to be involved in the catalytic 
reaction: K46, D130, K170, and E203 [16]. Alanine substitutions 
in the catalytic tetrad (K46, D130, K170, or E203) almost 
completely block 2'-0-MTase activity without jeopardizing 
binding to nsplO (Table 2, and [16]). Several SAM-binding 
residues (N43, G73, D99 and Y132, Fig. S2A) were substituted by 
alanine. Although they conserve their specific nsplO binding 
properties, indicating that they are correctly folded, they all show a 
drastically reduced MTase activity (Table 2), validating the 
structural description of the nsplO/nspl6/SAH ternary complex. 

A Mg 2+ Cation is Present in nspl 6, Outside the Active Site 

We recently reported[16] that nspl0/nspl6 MTase activity 
requires Mg 2+ . Although the crystallization buffer contains Mg 2+ , 
we were unable to locate any such cation in the nspl6 active site. 
In enzyme activity assays, the Mg 2+ ion can be substituted by 
Mn i+ or Ca 2+ ’ but not Zn 2+ (data not shown, see also [16]). A peak 
of electron density presumably corresponding to Mg 2+ is localized 
onto nspl6, distant from the SAH-binding cavity. The Mg 2+ 
coordination mode is through six first-shell water molecules in an 
octahedral geometry (Fig. 4). Binding via water molecules, involves 


T58 and SI88 side chain hydroxyls and the main chain carbonyl 
of E276. Since there are no carboxylic acids involved in binding 
this cation, it was suspected that its presence resulted from the 
crystallization procedure [42], with no biological relevance. 
However, the T58A, T58N, T58E and S188A substitutions show 
43, 70, 99 and 72% loss of activity, respectively (Table 2), with no 
significant effect on the stability of the nspl0/nspl6 complex 
except for T58E whose association was 54 % that of wild-type. 
These residues are located on three distinct structural elements at 
the C-terminus of helix Z (T58), the N-terminus of (36 (SI88), and 
in the central part of helix A3 (E276), respectively (Fig. 4 and Fig. 
S3). The cation may thus hold these elements together. 

The Putative RNA Binding Site: Mutagenesis and Effect 
on Activity 

The nspl0/nspl6 complex absolutely requires an N7-methyl 
guanine capped RNA substrate to exhibit MTase activity[16]. The 
structural basis for the preferential binding to methylated N7- 
guanine versus non-methylated caps has been elucidated in four 
cases, those of VP39[27], eIF4E[43], CBC[44], and PB2[45] 
proteins (PDB codes 1AV6, 1EJ1, 1H2T, and 2VQZ, respectively) 
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Figure 2. Topology diagram of MTase secondary structure elements. A) The consensus topology diagram is shown with a two domain 
organization for methyl acceptor substrate and SAM, as defined[25] with the catechol O-MTase and its canonical seven strand beta sheet. B) Topology 
diagram of nsp!6 colored according to rainbow colors from N-to C-terminus as in Fig. IB. The missing helix B and helix C are indicated by dashed 
ovals. The approximate general binding site of nspIO is shown (yellow oval). 
doi:10.1371/journal.ppat.1002059.g002 


bound to cap analogues or capped RNAs. In these cases, the 
methylated base specificity is achieved through increased 
binding energy resulting from the stacking of the N7-methyl 
guanine between parallel aromatic residues of the cap binding 
protein. The presence of the methyl group greatly enhances 7t-7l 
stacking, providing a dominant effect over unmethylated 
guanine [46]. Despite numerous attempts, cap analogues 
(m7GpppA, GpppA, m7GpppG, GpppG) and short capped 
RNA substrates (m7GpppA(C) n ) could neither be co-crystallized 
with nspl0/nspl6 nor soaked and bound onto preformed 
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nspl0/nspl6 crystals. However, the atomic coordinates of the 
N7-methyl guanine RNA oligomer in complex with VP39[47] 
provided data from which a model of RNA binding to the nsp 16 
protein was derived. SAM molecules identified in both 
structures were superimposed, and the VP39-bound RNA was 
positioned onto the nsp 16 structure. After minimal manual 
adjustments not exceeding 5 A, the VP39 RNA was a 
reasonably good fit into an nsp 16 hydrophobic groove radiating 
from the catalytic site (Fig. 5), establishing very few contacts 
with nsp 10. We note that the protein side diametrically opposite 
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Figure 3. The SAM binding site of nsp16. A) SAH modeled in a simulated annealing Omit map contoured at la. Carbons, oxygens, sulfur, and 
nitrogens are in grey, red, yellow, and blue sticks, respectively. Water molecules are shown as red spheres. B) Sinefungin bound in the SAM/SAH 
binding site, with main catalytic residues and the water molecule indicative of a catalytic mechanism. Colors as in A). C) SAH modeled in a 2Fo-Fc map 
(green) of Sinefungin contoured at la with a Fo-Fc difference map (red) contoured at 3a. A peak of negative density appears clearly on the sulfur 
atom of the SAH molecule showing that SAH was indeed replaced by Sinefungin. Colors as in A). 
doi:10.1371/journal.ppat.1002059.g003 


to the proposed hydrophobic RNA binding groove is highly 
positively charged (not shown), an observation that may account 
for the difficulty of achieving experimental RNA binding in the 
proposed RNA binding site. In the absence of robust data to 
guide docking of the guanine cap, the m7Gpp cap structure was 
not positioned in the structure but two possible N7-methylated 
cap guanine binding areas are indicated by arrows (Fig. 5). The 
first transcribed nucleotide together with its ribose receiving the 
methyl group fit well in the active site (Fig. 5, panel B) as 
predicted in the proposed mechanism. The same holds for the 
immediately preceding three nucleotides. The base of the first 
transcribed nucleotide may be held by contact with PI34 and 
Y132, bending the extending RNA cap structure. Accordingly, 
the substitution of Y132 greatly depresses MTase activity 
(Table 2). We also note that Y132 is located in the vicinity of 
a highly mobile loop (residues 135-138) not always visible in our 
crystal structures suggesting that this loop may move in order to 
wrap the triphosphate moiety of the RNA cap and/or the RNA 
cap itself. The solvent exposed side chain of Y30 may also 
participate in RNA binding. In the model, the highly mobile 
side chain of Y30 was flipped out in an alternative conformation 
in order to open the groove. In that position, Y30 should 
specifically contact the third transcribed nucleotide. Our 
mutagenesis data confirms the importance of Y30 since its 
replacement with either Ala or Phe severely impairs MTase 
activity without affecting the interaction with nsplO (Table 2). 

Interface of the Heterodimer 

All nsplO secondary structure elements but helices 2 and 5 
contact nspl6 (Fig. 1). The nsplO contact points can be viewed as 
5 small patches A to E (residues 40M7, 57-59, 69-72, 77-80, and 
93-96, respectively, Fig. 6A). In turn, these five patches contact 
most of the nsp 16 SAM-binding structural elements in 4 areas, I to 
IV (Fig. 6B, Fig. S3)), mainly involving (32, (33, ocA, aZ, and for 
area IV, the loop connecting helices A2 and A3 at the C-terminus 
(Figs. 1 and 6). In total, the interface of the heterodimer involves 
53 residues, 23 and 30 from nsplO and nsp 16, respectively. In 


nsplO, a single residue (AsnlO, at the edge of the interaction 
surface) is not conserved out of 23 (4.3 %), whereas in nsp 16 there 
are 8 non-conserved residues out of 30 (26.7 %) (Fig. SI). The 
interface has a buried surface area of 1820 A 2 , with nsplO 
contributing to 930 A 2 and nsp 16 to 890 A 2 . 

Four nsplO patches included in the 5 interaction patches 
identified here were recently mapped using reverse yeast two- 
hybrid methods coupled to bioluminescence resonance energy 
transfer and in vitro pull-down assays (see[18] and below). To probe 
the observed ciystal structure of the interface further, we 
engineered 5 new nsplO alanine mutants (N40A, L45A, T58A, 
G69A, and H80A, see Table SI) sitting in patches A, B, C and D. 
Whereas T58A, G69A and H80A showed limited effect on nsp 16 
binding, N40A reduced it to 64% of wild type affinity, and L45A 
almost abrogated it. The crystal structure indicates that the nsplO 
interface proposed by Lugari et o/.[18], is a correct and 
conservative estimation, as the interface also includes L45 
belonging to patch A. We also confirm the positive co-relation 
of the detected nspl0/nspl6 interaction with MTase activity. In 
no instance can nsp 16 be active in the absence of nsp 10/16 
complex formation. 

The Y96 position is of particular interest. Alanine substitution 
(Y96A) abrogates interaction whereas a phenylalanine (Y96F) 
increases both interaction and MTase activity[18]. In order to 
understand how residue 96 plays such a pivotal role, we 
determined at 2.0 A resolution the crystal structure of this 
nsplO(Y96F)/nspl6 complex (Table 1). Strikingly, the absence of 
the hydroxyl group does not alter the topology of the interface. 
Wild-type and Y96F residues superimpose without significant 
difference at all atomic positions (not shown). Either Y96 or F96 is 
in direct contact with nsp 16 helix aZ, which carries the catalytic 
residue K46. Detailed surface analysis using PISA indicates that 
the position of K46 in Y96F nsplO is identical to that of K46 in 
wild-type nsp 10, ruling out a better alignment of catalytic residues 
of the Y96F mutant. The nsplO(Y96F)/nspl6 differs from wild- 
type nspl0/nspl6 in the SAH binding site, though. When 
compared to that of wild-type, the SAH occupancy is much lower 



PLoS Pathogens | www.plospathogens.org 


6 


May 2011 | Volume 7 | Issue 5 | el 002059 




SARS-CoV nsp10/16 MTase Complex Structure 


Table 2. Mutational analysis, complex formation, and enzyme activity of the nsp10/nsp16 complex. 


Number 

Function* 

Mutant 

% associated** 

% MTase act. *** 

1 = Wild-Type 

Wild-TypE 

Wild-Type 

100 

100 

2 

Catalytic 

K46A 

97 

1 

3 

Catalytic 

D130A 

104 

2 

4 

Catalytic 

K170A 

95 

0 

5 

Catalytic 

E203A 

93 

0 

6 

SAM-BS 

N43A 

144 

11 

7 

SAM-BS 

G73A 

119 

18 

8 

SAM-BS 

D99A 

57 

0 

9 

Mg 2+ -BS 

T58A 

123 

57 

10 

Mg 2+ -BS 

T58N 

110 

30 

11 

Mg 2+ -BS 

T58E 

54 

1 

12 

Mg 2+ -BS 

SI 88 A 

102 

28 

13 

RNA/SAM-BS 

Y132A 

86 

5 

14 

RNA/SAM-BS 

Y132T 

95 

5 

15 

RNA/SAM-BS 

Y132F 

123 

9 

16 

RNA/SAM-BS 

Y132H 

106 

0 

17 

RNA-BS 

Y30A 

89 

1 

18 

RNA-BS 

Y30F 

113 

6 

Interface mutants 

Number 

Location* 

Mutant 

% associated** 

% MTase act. *** 

19 

Patch 1 

I40A 

84 

8 

20 

Patch 1 

M41A 

122 

4 

21 

Patch 1 

V44A 

113 

1 

22 

Patch 1 

T48A 

123 

20 

23 

Patch II 

V78A 

3 

0 

24 

Patch II 

R86A 

61 

0 

25 

Patch II 

Q87A 

114 

62 

26 

Patch III 

V104G 

6 

4 

27 

Patch III 

D106A 

101 

38 

28 

Patch IV 

L244A 

30 

0 

29 

Patch IV 

M247A 

4 

0 

*as inferred from the crystal structure. 


**nsp10/nsp16 complex formation as determined using pull-down experiments, SDS-PAGE analysis, and quantitation (see Methods). 
***The % of MTase activity was determined using filter binding assays (see Methods) relative to wild-type. 

SAM-BS: S-adenosylmethionine binding site. 

Mg 2+ ~BS: Magnesium binding site. 

RNA-BS: RNA binding site. 

doi:10.1371/journal.ppat.1002059.t002 


(~0.3 versus ~1), leading to poor density definition. SAH is 
known to be a fairly good inhibitor of the methylation reaction. 
Therefore, a lower binding affinity might translate into less end- 
product inhibition, and account for the observed increased 
activity. We measured the affinity of nsplb for SAH using 
fluorescence spectroscopy, but no significant differences were 
found (not shown). Likewise, the MTase inhibition pattern by 
SAH was identical for wild-type and nsplO(Y96F)/nspl6 (not 
shown). We therefore infer that the previously observed ~ 10-fold 
increased stability of the heterodimer [18] may be responsible for 
the increased activity. A more hydrophobic character of the 
interaction may appear upon the loss of the tyrosine hydroxyl 
which, in the wild-type protein, was not engaged in any polar 
contact. We thus attribute the increased activity of the 
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nsplO(Y96F)/nspl6 complex relative to wild-type to a stronger 
equilibrium association of nsplO(Y96F) with nspl6 than that of 
wild-type nsplO with nspl6. 

Mutation analysis was also conducted on nspl6 residues 
presumably involved in the interface and interfacial activation 
(Tables 2 and S2). Several mutants (V78A, V104A, L244A, 
M247A) in patches II, III and IV completely disrupt the nsplO/ 
nsp 16 complex and annihilate nsp 16 MTase activity. Interestingly, 
we also identified nsp 16 mutants still interacting with nsplO, but 
with a strongly reduced 2'-0-MTase activity (I40A, M41A, V44A, 
T48A, Q87A, D106A) suggesting, that these mutations in the 
nspl0/nspl6 interface may alter the fine positioning of catalytic 
residues without any significant effect on nsplO binding. 
Accordingly, most of these mutants are localized in aZ helix of 
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Figure 4. Position and coordination of a metal ion. Left, the putative Mg 2+ ion (blue sphere) is shown solvated in its first atomic shell by six 
water molecules (red spheres). The corresponding 2Fo-Fc electron density map (contoured la) with a cross shape is shown in purple. Residues of 
Nsp16 involved in coordination of Mg 2+ ion via water molecules are labeled. Right, a global view of the bound metal ion in its electron density on the 
opposite side from the SAFI molecule (red sticks). Nspl 6 is shown as green ribbons, nspl0 as yellow ribbons with its two bound Zn 2+ as grey spheres 
as in Fig. 1. 

doi:10.1371/journal.ppat.l 002059,g004 


patch I which contains the K46 catalytic residue. On the other 
hand, patch II and III mutants tend to have more mitigated 
phenotypes, yielding to full-blown interaction with only about half 
of the expected activity. Finally, patch IV mutants were totally 
inactive. Using all mutants reported in Table 2, a plot (Fig. 6C) of 
interaction versus activity shows that the nsp 10/nsp 16 interaction 
is strictly required to obtain significant nsp 16 MTase activity. 

Discussion 

Mechanism of 2'-0-Methyltransfer 

The SARS-CoV RNA cap 2’-0-MTase is a heterodimer 
comprising SARS-CoV nsp 10 and nsp 16. When bound to nsp 10, 
nspl6 is active as a type-0 RNA cap-dependent 2'-0-MTase, ie., 
active only when the cap guanine is methylated at its N7 
position[16]. The nspl0/nspl6 crystal structure shows that nspl6 
adopts a typical fold of the S-adenosylmethionine-dependent 
methyltransferase family as defined initially for the catechol O- 
MTase[25]. A good alignment (170°) is found between the SAH 
sulfur atom, a water molecule present in both SAFI- and sinefungin- 
bound nsp 16 structures, and the K46 8-amino group (Fig. 3A and 
B). This geometry provides interesting hints for a catalytic 


mechanism, as the positions of the catalytic residues (K46, D130, 
K170, E203) match spatially those of the vaccinia virus VP39 2’-0- 
MTase[47], At the initial stage of the reaction the 2 , -hydroxyl of the 
capped RNA substrate would occupy the position of die water 
molecule. In turn, E203 and K170 decrease the pKa of the K46 s- 
amino group that becomes a deprotonated general base (-NH 2 ) able 
to activate the RNA 2 , -hydroxyl at neutral pH. In VP39, K175 has 
been identified as the general base catalyst[47] with a pKa 
depressed by ~ 2 pH units by the neighbouring D138 and R209 
residues [48], These findings indeed suggest a related mechanism: 
once K46 has activated the 2'-hydroxyl group, the 2'-oxygen would 
produce an in line attack through a SN2-like mechanism onto the 
electrophilic SAM methyl group. The methyl group would pass 
through a pentavalent intermediate with the 2’-0 and sulfur at 
apical positions. D130 is positioned to stabilize the transient positive 
charge on the donated methyl atom of SAM before the sulfur 
recovers a neutral electric charge during SAH generation (Fig. 3B). 

Role of Mg 2+ on MTase Activity 

Unlike most SAM-dependent MTases, the SARS-CoV nsplO/ 
nsp 16 enzyme requires a divalent cation, either magnesium, 
manganese or calcium[25]. We have found that this cation does 
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Figure 5. Stick model of RNA bound to the nsp16 RNA binding groove and Sinefungin in the methyltransferase active site. A) In this 
representation, Sinefungin was preferred over SAH because one of its NH 2 groups approximates the direction of a transferred CH 3 group from the 
SAM substrate. Carbon is white, oxygen is red, nitrogen is blue and phosphorous is orange. Nsp16 and nspIO are rendered as a solvent-accessible 
surface colored grey and wheat respectively. The Sinefungin molecule defines the methyltransferase active site. Missing residues in the 135-137 loop 
(see text) are indicated by a shaded blue dotted box. Position of Y30 and Y132 are indicated. Y30 generated poor electron density (see « Methods ») 
and its aromatic ring position has been manually adjusted before generation of this image. B) Close caption of the methyltransferase active site 
showing distance between the NH2 of Sinefungin to the 2'-0 of the ribose of the first base, thus mimicking the position of the methyl of the S- 
adenosylmethionine. 
doi:10.1371/journal.ppat.1002059.g005 


not reside in the active site. Instead, the cation is coordinated 
through water molecules by three residues located on three distinct 
structural elements. It is thus possible that one divalent cation, 
presumably Mg~ + , present in the host cell at miUimolar levels, 
plays a structural role in holding these three nspl6 structural 
elements together and so regulate the enzyme activity. It is 
intriguing that T58A is more active than T58N or T58E that can 
still bind the water that chelates to the metal. Alternatively, it is 
possible that divalent cations such as Mg z+ or Ca z+ act as a 
phosphodiester charge shield to allow RNA binding in the 
hydrophobic binding groove [49]. 

MTase Activity is Regulated via Protein-Protein 
Interactions 

The main regulation mechanism of nsp 16 is through its physical 
association with nsp 10. Nsp 16 is unstable in solution, and nspIO 
acts as a scaffold for nsp 16, yielding a stable dimer active as an 
RNA cap-dependent (nucleoside-2'-0)-MTase. The complex is 
assembled through a ~890 A z contact surface in nsp 16, an area 
typically in the intermediate zone differentiating strongly from 
weakly associated dimers [50], This finding is consistent with a Kd 
estimated at ~0.8 |J,M[18] that qualifies the nspl0/nspl6 complex 
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as a rather weak heterodimer. The nspIO interaction surface 
identified in the crystal structure was confirmed by site-directed 
mutagenesis and overlaps that previously identified by indirect 
methods[18]. Remarkably, the nspIO surface in the nspl0/nspl6 
complex is essentially identical to that of uncomplexed nspIO 
crystallized alone by others [20,21] (Fig. S4). It is therefore 
reasonable to see this heterodimer as a non-permanent species 
which would tolerate nsp 10 or nsp 16 engaging in interactions with 
other partners. This notion is actually in line with the involvement 
of nsp 10 in a network of protein-protein interactions that we and 
others have proposed[17,19]. Donaldson et a/.[23] have engi¬ 
neered mutations in nspIO using reverse genetics. Out of eight 
mutations that turned out to be in the nspl0/nspl6 interface (this 
work), five, two and one rendered lethal, debilitated, and viable 
phenotypes, respectively[23]. Interestingly, the nsplO(Q65E) 
mutant providing a temperature-sensitive phenotype[22,23] does 
not map in the nspl0/nspl6 interface, confirming that nspIO has 
a pleiotropic role. 

Our mutagenesis analysis shows that the formation of an nspIO/ 
nsp 16 complex is a pre-requisite for MTase activity (Fig. 6C) 
indicating that physical association of nsp 10 and nsp 16 is essential 
to activate nspl6 2'-0-MTase activity and foster efficient virus 
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Figure 6. Detailed definition of nsp16 and nspIO elements involved in the interface. A) nsp16 (left) and nspIO (right) are rendered as a 
solvent-accessible surface colored green and yellow, respectively. Residues involved in the interface are rendered as sticks with the same color code, 
in transparency into respective proteins. B) Separate representation of the nsp16 (left) and nspIO (right) interface. Nspl6 residues defining patch I, II, 
III, and IV are rendered as sticks as above in transparency into nsp16. NspIO residues defining patch A, B, C, D, and E are rendered as sticks as above in 
transparency into nspIO. C) A plot of MTase activity as a function of nsp10/nsp16 complex formation using the data reported in Table 2. 
doi:10.1371 /journal.ppat.1002059,g006 


replication. We note that most interface mutants exhibit a severe 
loss of their 2'-0-MTase activity, whereas the apparent associa¬ 
tion affinity is often only modestly affected. That minor changes in 
the interface translate into potent effects is also dramatically 
illustrated by the Y96F mutation, where the loss of a single 
hydroxyl provokes a significant change in affinity[18]. Remark¬ 
ably, it is not the most active complex that was selected in nature, 
since the nsplO(Y96F)/nspl6 complex is both more stable and 
more active than the wild-type heterodimer (this work and[18]). 
This is yet another observation hinting at the involvement of nsp 10 
in protein-protein interaction networks including other partners 
than nsp 16, such as nsp5 and nsp 14[17,19]. In most other 
coronaviruses, the nspIO residue at position 96 is a phenylalanine. 
It would be interesting to determine whether this polymorphism is 
relevant to the SARS-CoV pathogenicity at any (direct or indirect) 
level, or if compensating polymorphisms in other coronaviral 
nspIO (or nsp 16) restore a weaker nspl0/nspl6 association 
equivalent to that of the SARS-CoV pair. Since a bona fide viral 
RNA cap is key in evading the host cell innate immunity[4,51], a 
minimal level of 2'-0-MTase activity would be expected to be 
critical to virus survival. 

MTase activation through dimerisation of two viral protein 
partners has already been reported in the case of the vaccinia virus 
D1/D12 N7-guanine MTase[28]. However, the activating D12 
subunit does not contact the D1 subunit through a homologous 
surface mainly defined by canonical ctA and aZ helices. Rather, 
the D1 /D12 activation surface would be located at a 90° clockwise 
rotation relative to the nspl0/nspl6 interface depicted in Fig. 1A. 
In the case of dengue virus, the bi-functional N7-guanine and 2'- 
O-MTase is part of the N-terminus of the dengue NS5 protein. 
Based on reverse genetic data and modeling[52], the MTase 
domain would be associated with the Pol domain through an 
interface topologically similar to that of nspl0/nspl6, i.e., 
involving mainly helices ocA, aZ and strands [12 and [13 as 
depicted in Fig. 1A. 

We have previously shown that the nspl0/nspl6 is only active 
as N7-guanine methylated capped RNA, implying that RNA cap 
methylation obeys to an ordered sequence of events where nsp 14- 
mediated N7-guanine methylation precedes nspl0/nspl6 RNA 
2'-0 methylation[16]. In the absence of data regarding the RNA 
substrate, we built a model of RNA binding based on that of the 
vaccinia virus VP39 ternary complex structure. Interestingly, our 
model proposes that the RNA interacts only with nsp 16 residues, 
in keeping with what was recently suggested based on RNA 
binding assays [14], Although the position of the cap structure on 
the nsp 16 surface remains to be determined, our model suggests a 
well-defined position for the ribose of the first transcribed 
nucleotide in the active site. In agreement with mutagenesis 
analysis, the model also suggests that the transcribed RNA 5'-end 
stacks between Y132 and Y30. Furthermore, this model is 
consistent with the observation that coronavirus MTase requires 
RNA substrates of at least 3 transcribed nucleotides in length [14], 
It is also worth to know that a comparison of nsp 16 and VP39 
electrostatic surfaces reveals that the putative RNA-binding groove 
of nsp 16 is mostly hydrophobic, whereas the VP39 RNA-binding 
groove is positively charged. This variation would imply a change 
in the nature of the RNA/protein interaction. 


Binding of Sinefungin to nspl6/nsp!0 Suggests 
Feasibility of a Drug Design Approach to Inhibit MTase 
Activity 

Viral MTases are increasingly evaluated as potential drug 
design targets[34,37,53]. We have crystallized the inhibitor 
Sinefungin with the nspl0/nspl6 complex. Sinefungin exhibits 
an IC 50 of 0.74 |xM, 16-fold lower than that of SAH as reported by 
Bouvet et al. [16] using purified nspl0/nspl6. Analysis of the 
structure suggests a likely mechanism of action that also accounts 
for the observed inhibitory effect of diis drug. We note that the 
adenine nucleobase does not fit snugly into its binding pocket, 
raising interest regarding structure-based drug design. Preliminary 
examination of eukaryotic non-viral MTase structures from main 
classes as defined in Martin and McMillan [25] indicates that the 
SAH adenine is bound tighter in any of the latter enzymes than in 
the nsp 16 SAM-binding site, indicating a possible breach to 
achieve anti-coronavirus selectivity with a small molecule inhibitor 
of nsp 16. 

In conclusion, the crystal structures presented here extend our 
general understanding of the mechanism and regulation of viral 
RNA cap MTases in (+)RNA viruses, and point to both the nspIO/ 
nsp 16 interface and the substrate binding sites as putative antiviral 
targets. 

Methods 

Crystallization and Data Collection 

Both nspIO and nsp 16 were expressed from the same dual 
expression vector pmCOX [16]. NspIO had a N-terminal strep-tag 
(WSHPQFEK), and nsp 16 a N-terminal hexa-histidine tag. The 
purification and crystallogenesis of the nspl0/nspl6 complex was 
performed as described in [40], Typical crystals of the wild-type 
nspl0/nspl6 appear in hanging drops after 24 h at 20°C in 0.1 M 
CHES pH 9, 1.52 M MgCE hexahydrate. Crystals (a = 68.53 A, 
b= 184.74 A, c= 129.01 A, C222i) contain one nspl0/nspl6 
complex per asymmetric unit, with a solvent content of 70 % and 
V m of 4.17 A 3 /Da. Crystals of nsplO(Y96F)/nspl6 were grown in 
67 mM CHES pH 8.5, 0.99 M MgCE hexahydrate, 33 mM Tris- 
HC1, 8.3 % PEG 8000. Both crystallization conditions yielded 
crystals diffracting to 1.9 A when exposed to synchroton radiation 
at the ID 14-1 beamline of the European Synchrotron Radiation 
Facility, Grenoble, France. Crystals were cryo-cooled in the same 
buffer supplemented with 15 % glycerol. Crystal soaking was 
performed in the same buffer supplemented with 5 mM SAH or 
Sinefungin during 24 h. 

Structure Determination and Refinement 

The position of the nspIO protein was unambiguously 
determined by molecular replacement using the program 
PHASER[54] with the nspIO protein (2FYG), as search probe[20]. 
Strong peaks in both the residual and anomalous Fourier maps 
confirmed the presence of two Zinc ions at the expected positions 
within the nsp 10 protein, thus giving confidence in the validity of 
the MR solution. Phases calculated from this partial model were 
combined with SAD phases from the Zn atoms using PHASER. 
To ameliorate the resulting low quality density map, phases were 
improved with PARROT [55]. An initial model, comprising both 



PLoS Pathogens | www.plospathogens.org 


11 


May 2011 | Volume 7 | Issue 5 | el 002059 



SARS-CoV nsp10/16 MTase Complex Structure 


nsplO and nspl6, was automatic built by successive use of 
BUCCANEER[56] and ARP/wARP[57], The resulting model 
was subject to several cycles of manual rebuilding using COOT 
[58] and refinement with REFMAC [59]. The protein structure 
model could be built, except the strep and hexahistidine tags. In 
nspl6, density was too weak for the mobile, solvent exposed nspl6 
loop 136-139, Y30 (see “Results”), and 2 and 6 residues in N- and 
C-terminus, respectively. Likewise, nsplO solvent exposed 9 and 8 
residues in N- and C-terminus were missing, respectively. 

Overall, the chain traces are unambiguous, with clear electron 
density including for a single SAH residue bound to the nspl6 
protein. Solvent accessible surfaces were calculated using program 
AREAIMOL [60] with a 1.7 A radius sphere as the probe 
(Table 1) and values rounded to the nearest 5 A 2 . Conformational 
differences were analyzed using the DynDom server (http://www. 
cmp.uea.ac.uk/dyndom/main.jsp). Figures were created using 
PYMOL (http://www.pymol.org). The coordinates of the wild- 
type/SAH, mutant, and wild-type/Sinefungin structures have 
been deposited at the Protein Data bank under PDB codes 2XYQ, 
2XYV, and 2XYR, respectively. 

RNA Cap Structure Modeling 

The modeling of the RNA cap structure in the nspl0/nspl6 
complex structure is derived from the analysis of the structure of 
the vaccinia virus methyltransferase VP39 crystallized in complex 
with a capped RNA and a S-Adenosylhomocysteine [47] (SAH) 
(pdb code: 1AV6). The two structures are manually aligned using 
COOT [58] based on the position of SAH binding sites, as well as 
SAH, and Sinefungin (SFG) molecules. The RNA binding site of 
VP39 is only partly overlapping that of nsp 16 whilst the shape of 
the cavity is similar; thus local adjustments necessary to 
accommodate the RNA molecule in its binding groove were done 
manually using COOT. The side chain of tyrosine 30 of 
nsp 16initially pointed to the putative RNA binding site, preventing 
any bona fide modeling. In order to fit the RNA molecule in the 
cavity, an alternative conformation was sought for this side chain. 
The second most common conformation for the tyrosine side 
chain was selected. Due to the biochemistry data and surface 
electrostatic analysis, it is not possible to describe with certainty the 
final position of the cap, thus the cap was removed and replaced 
by arrows symbolizing possible positions. No other modification 
was performed on the RNA, the Sinefungin molecule or the nsp 16 
structure. 

Plasmids 

The SARS-CoV nsplO and nsp 16-coding sequences were 
amplified by RT-PCR from the genome of SARS-CoV Frank¬ 
furt-1 (accession number AY291315) as previously described [16]. 
The nsplO and nspl6 genes (encoding residues 4231-4369, 5903- 
6429, and 6776-7073 of replicase pplab) were cloned into a 
Gateway modified dual-promotor expression plasmid and in the 
gateway pDest 14 expression vector. In this backbone, SARS CoV 
nsplO can be expressed under a tet promoter and encodes a 
protein in fusion with a N-terminal strep tag, whereas nsp 16 is 
expressed under a T7 promoter and encodes a protein in fusion 
with a N-terminal hexahistidine tag. The mutants were generated 
by PCR using the Quickchange site-directed mutagenesis kit 
(Stratagene), according to the manufacturer’s instructions. 

Reagents 

AdoMet and cap analogs GpppA and 7Mc GpppA were 
purchased from New England BioLabs, the [’H]-AdoMet was 
purchased from Perkin Elmer and Sinefungin (adenosylornithine) 
from Sigma-Aldrich. 


Expression and Purification of SARS-CoV nsplO, nsp!6 
and nsp!0/nsp!6 Complex 

E. coli C41 (DE3) cells (Avidis SA, France), containing the pLysS 
plasmid (Novagen), were transformed with nsplO or nsp 16 cloned 
in pDest 14, or nspl0/nspl6 cloned in pmCox, and grown in 2YT 
medium supplemented with appropriate antibiotics. The expres¬ 
sion of strep-tagged nsplO or 6His-tagged nsp 16 mutants was 
induced (DO 600 = 0.6) by adding 50 pM IPTG, and the expression 
of the nsp 10/nsp 16 complex by adding 50 pM IPTG and 200 pg/ 
L of anhydrotetracycline. After an incubation for a 16 h at 24°C, 
the cell were pellets, frozen and resuspended in lysis buffer 
(50 mM HEPES, pH 7.5, 300 mM NaCl, 5 mM MgS0 4 , 5 mM 
fi-mercaptoethanol (only for nsplO) supplemented with 1 mM 
PMSF, 40 mM imidazole, 10 pg/ml DNase I, and 0.5% Triton 
X-100. After sonication and clarification, proteins were purified 
either by IMAC (HisPurTM Cobalt Resin; Thermo Scientific) 
chromatography[16] (nsplO mutants and nspl6 mutants), and the 
nsp 10/nsp 16 complex was purified by using Strep-Tactin 
sepharose (IBA Biotagnology) as previously described [16], All 
purified proteins were analyzed by SDS-PAGE. The binding of 
wild-type nsplO to mutant nsp 16, and that of mutant nsplO to 
wild-type nsp 16 was quantified using ImageJ as described [16]. 

Radioactive Methyltransferase Assay 

MTase activity assays were performed in 40 mM Tris-HCl, 
pH 8.0, 5 mM DTT, 1 mM MgCl 2 , 2 pM 7Mc GpppAC 5 or 
GpppAC 5 , 10 pM AdoMet, and 0.03 pCi/pl [ ,5 H] AdoMet (GE 
Healthcare). Short capped RNAs ( 7M ‘‘GpppAC 5 , GpppAC 5 , were 
synthesized in vitro using bacteriophage T7 DNA primase and were 
purified by high-performance liquid chromatography (HPLG) as 
previously described[61]. In the standard assay, nsplO and nsp 16 
were added at final concentrations of 600 nM, and 200 nM, 
respectively, and the amount of ’H-GH :5 transferred onto 
7Ml 'GpppAC 5 substrates was determined by filter binding assay 
as previously described [16], 

Supporting Information 

Figure SI Sequence alignment and amino acid conservation in 
coronavirus nsp 16. The alignment of coronavirus nsp 16 sequences 
was generated with Muscle program (http://www.ebi.ac.uk/Tools/ 
msa/muscle/), and the resulting alignment converted using the 
ESPript program, (http://espript.ibcp.fr/ESPript/cgi-bin/ESPript. 
cgi). Residues that are conserved in all or >70% sequences are 
boxed in red and yellow, respectively. National Center for 
Biotechnology Information (NCBI) accession numbers for replicase 
polyprotein sequences that include nsp 16 are as follows: Severe 
acute respiratory syndrome virus SARS-CoV Frankfurt isolate 
(SARF), Severe acute respiratory syndrome virus SARS-CoV Tor2 
isolate (SART), NP_828873.2; Turkey coronavirus (Turk), 
YP_001941189; Infectious Bronchitis Virus (IBV), NP_066134; 
Feline Coronavirus (Feli), YP_239426;Porcine Transmissible Gas¬ 
troenteritis Coronavirus (PTGC), P18457; Transmissible Gastro¬ 
enteritis coronavirus (TGEV), NP_840013; Porcine epidemic 
diarrhea virus CV777 (PEDV), NP_839969; Bat coronavirus 
512/2005 (BATC), YP_001351683; Human coronavirus NL63 
(NL63), AAS58176; Human coronavirus 229E (229E), NP_073549; 
Murine hepatitis virus strain JHM (MuiJ), YP_209243; Mouse 
Hepatitis virus strain A59 (MA59), NP_740613; Human coronavi¬ 
rus HKU1 genotype A (HKU1), YP_460023; Human enteric 
coronavirus 4408 (4408), ACJ35483; Bovine coronavirus (BCoV), 
NP_742142; Human coronavirus OC43 (OC43), AAT84359. 
Numbering was made using SARF as a reference. 

(TIF) 
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Figure S2 Binding determinants in the SAM/SAH/Sinfungin 
binding site. A) LIGPLOT diagram (www.ebi.ac.uk/thornton-srv/ 
software/LIGPLOT/) of the SAH ligand molecule interacting 
with the nspl6 binding site. Ligand bonds are in purple, neighbor 
residue (non-ligand) bonds are in light brown, hydrogen bonds are 
green dashed lines. Ligand atoms surrounded by a yellow circle 
are highly accessible. Non-ligand residues in hydrophobic contacts 
with the ligand are presented by red semi-circles with radiating 
spokes. B) Superimposition of SAH and Sinefungin molecules. The 
nspl6 residues are green sticks, the SAH is represented in red 
sticks, and Sinefungin and water molecules are in sticks and 
spheres as in Fig. 3, respectively. 

(TIF) 

Figure S3 Primary and secondary structure elements of nspl6. 
Helices and sheets are colored according to rainbow colors from 
N- to C-terminus as in Fig. 1. Patches I to IV are boxed in grey, 
labeled above the sequences, and correspond to residues 
represented in Figs. 6A and B. 

(TIF) 

Figure S4 NsplO interface comparison between the nsplO/ 
nspl6 complex and nsplO homo multimer. A) For the corre¬ 
sponding nsplO/nsplb interface nsplO (cyan) interacts with 2 
monomers (green and orange). Left: Cartoon representation of 3 
nsplO monomer part of the dodecamer structure (PDB: 2G9T). 
Right: Cartoon representation of the nsp 10(cyan)-nsp 16(dark red) 
complex. NsplO molecules in cyan are in the same orientation. B) 
NsplO sequence is presented with above the corresponding 
secondary structure elements. Below dark red colored dots indicate 
residues involved in the interaction between nsplO and nsp 16 
while green and orange dots indicate residues involved in the 
homomultimer in the dodecamer complex. 

(TIF) 

References 

1. Gu M, Lima CD (2005) Processing the message: structural insights into capping 
and decapping mRNA. Curr Opin Struct Biol 15: 99-106. 

2. Shuman S (2001) Structure, mechanism, and evolution of the mRNA capping 
apparatus. Prog Nucleic Acid Res Mol Biol 66: 1—40. 

3. Shuman S (2002) What messenger RNA capping tells us about eukaryotic 
evolution. Nat Rev Mol Cell Biol 3: 619-625. 

4. Yoneyama M, Fujita T (2010) Recognition of viral nucleic acids in innate 
immunity. Rev Med Virol 20: 4-22. 

5. Shuman S (2001) The mRNA capping apparatus as drug target and guide to 
eukaryotic phylogeny. Cold Spring Harb Symp Quant Biol 66: 301—312. 

6. Ahola T, Kaariainen L (1995) Reaction in alphavirus mRNA capping: 
formation of a covalent complex of nonstructural protein nsPl with 7-methyl- 
GMP. Proc Natl Acad Sci U S A 92: 507-511. 

7. Ogino T, Banerjee AK (2007) Unconventional mechanism of mRNA capping 
by the RNA-dependent RNA polymerase of vesicular stomatitis virus. Mol Cell 
25: 85-97. 

8. Plotch SJ, Bouloy M, Ulmanen I, Krug RM (1981) A unique cap(m7GpppXm)- 
dependent influenza virion endonuclease cleaves capped RNAs to generate the 
primers that initiate viral RNA transcription. Cell 23: 847—858. 

9. Rota PA, Oberste MS, Monroe SS, Nix WA, Campagnoli R, et al. (2003) 
Characterization of a novel coronavirus associated with severe acute respiratory 
syndrome. Science 300: 1394—1399. 

10. Gorbalenya AE, Enjuanes L, Ziebuhr J, Snijder EJ (2006) Nidovirales: evolving 
the largest RNA virus genome. Virus Res 117: 17—37. 

11. Lai MM, Patton CD, Stohlman SA (1982) Further characterization of mRNA’s of 
mouse hepatitis virus: presence of common 5'-end nucleotides. J Virol 41: 557-565. 

12. Lai MM, Stohlman SA (1981) Comparative analysis of RNA genomes of mouse 
hepatitis viruses. J Virol 38: 661-670. 

13. van Vliet AL, Smits SL, Rottier PJ, de Groot RJ (2002) Discontinuous and non- 
discontinuous subgenomic RNA transcription in a nidovirus. EMBO J 21: 
6571-6580. 

14. Decroly E, Imbert I, Coutard B, Bouvet M, Selisko B, et al. (2008) Coronavirus 
nonstructural protein 16 is a cap-0 binding enzyme possessing (nucleoside-2'O)- 
methyltransferase activity. J Virol 82: 8071-8084. 

15. Chen Y, Cai H, Pan J, Xiang N, Tien P, et al. (2009) Functional screen reveals 
SARS coronavirus nonstructural protein nsp 14 as a novel cap N7-methyltrans- 
ferase. Proc Natl Acad Sci U S A 106: 3484-3489. 


Table SI Effect of alanine mutations in nsplO interface residues. 
Residues were identified using PISA (http://www.ebi.ac.uk/ 
msd-srv/prot_int/pi_link.html). Bold on grey background: strictly 
conserved residues amongst coronaviruses; Bold on white back¬ 
ground: conserved a.a. (>70%) amongst coronaviruses (see Fig. SI). 
The % of Bioluminescence Resonance energy Transfer (BRET) 
signal was previously reported [18]. The interaction of each nsp 10 or 
nsp 16 mutant was determined using strep-tactin pull-down 
experiments of strep-tagged nsplO co-expressed with nsp 16 
followed by SDS-PAGE analysis, and quantitation (see Methods). 
The interaction of wild-type nsplO with wild-type nsp 16 was 
normalized to 100%. The % of MTase activity was determined 
using filter binding assays (see Methods) relative to wild-type. 
(DOC) 

Table S2 Effect of alanine mutations in nsp 16 interface 
residues. Residues were identified using PISA (http://www.ebi. 
ac.uk/msd-srv/prot_int/pi_link.html). Same legend as Table SI 
except that no BRET experiments were performed. 

(DOC) 

Acknowledgments 

We thank Barbara Selisko for helpful discussion, comments and support 
along this work, and Joelle Boretto-Soler for help in mutagenesis and 
protein purification. We are indebted to Karen Dalle, Violaine Lantez for 
excellent technical assistance and General Hugues Toulou for his support. 

Author Contributions 

Conceived and designed the experiments: E. Decroly, B. Canard. 
Performed the experiments: E. Decroly, C. Debarnot, F. Ferron, M. 
Bouvet, B. Coutard, I. Imbert, L. Gluais, N. Papageorgiou, B. Canard. 
Analyzed the data: E. Decroly, F. Ferron, I. Imbert, N. Papageorgiou, A. 
SharfF, G. Bricogne, M. Ortiz-Lombardia, J. Lescar, B. Canard. Wrote the 
paper: E. Decroly, B. Canard. 


16. Bouvet M, Debarnot C, Imbert I, Selisko B, Snijder EJ, et al. (2010) In vitro 
reconstitution of SARS-coronavirus mRNA cap methylation. PLoS Pathog 6: 
el000863. 

17. Imbert I, Snijder EJ, Dimitrova M, Guillemot JC, et al. (2008) The SARS- 
Coronavirus PLnc domain of nsp3 as a replication/transcription scaffolding 
protein. Virus Res 133: 136-148. 

18. Lugari A, Betzi S, Decroly E, Bonnaud E, Hermant A, et al. (2010) Molecular 
mapping of the RNA Cap 2'-0-methyltransferase activation interface between 
SARS coronavirus nsplO and nsp 16. J Biol Chem 285: 33230-33241. 

19. Pan J, Peng X, Gao Y, Li Z, Lu X, et al. (2008) Genome-wide analysis of 
protein-protein interactions and involvement of viral proteins in SARS-CoV 
replication. PLoS One 3: e3299. 

20. Joseph JS, Saikatendu KS, Subramanian V, Neuman BW, Brooun A, et al. 
(2006) Crystal structure of nonstructural protein 10 from the severe acute 
respiratory syndrome coronavirus reveals a novel fold with two zinc-binding 
motifs. J Virol 80: 7894-7901. 

21. Su D, Lou Z, Sun F, Zhai Y, Yang H, et al. (2006) Dodecamer structure of 
severe acute respiratory syndrome coronavirus nonstructural protein nsplO. 
J Virol 80: 7902-7908. 

22. Sawicki SG, Sawicki DL, Younker D, Meyer Y, Thiel V, et al. (2005) Functional 
and genetic analysis of coronavirus replicase-transcriptase proteins. PLoS Pathog 
1: e39. 

23. Donaldson EF, Sims AC, Graham RL, Denison MR, Baric RS (2007) Murine 
hepatitis virus replicase protein nsplO is a critical regulator of viral RNA 
synthesis. J Virol 81: 6356-6368. 

24. Donaldson EF, Graham RL, Sims AC, Denison MR, Baric RS (2007) Analysis 
of murine hepatitis virus strain A59 temperature-sensitive mutant TS-LA6 
suggests that nsplO plays a critical role in polyprotein processing. J Virol 81: 
7086-7098. 

25. Martin JL, McMillan FM (2002) SAM (dependent) I AM: the S-adenosylme- 
thionine-dependent methyltransferase fold. Curr Opin Struct Biol 12: 783—793. 

26. Vidgren J, Svensson LA, Liljas A (1994) Crystal structure of catechol O- 
methyltransferase. Nature 368: 354-358. 

27. Hodel AE, Gershon PD, Shi X, Quiocho FA (1996) The 1.85 A structure of 
vaccinia protein VP39: a bifunctional enzyme that participates in the 
modification of both mRNA ends. Cell 85: 247-256. 



PLoS Pathogens | www.plospathogens.org 


13 


May 2011 | Volume 7 | Issue 5 | el 002059 



SARS-CoV nsp10/16 MTase Complex Structure 


28. De la Pena M, Kyrieleis OJ, Cusack S (2007) Structural insights into the 
mechanism and evolution of the vaccinia virus mRNA cap N7 methyl- 
transferase. EMBO J 26: 4913-4925. 

29. Reinisch KM, Nibert ML, Harrison SC (2000) Structure of the reovirus core at 
3.6 A resolution. Nature 404: 960-967. 

30. Sutton G, Grimes JM, Stuart DI, Roy P (2007) Bluetongue virus VP4 is an 
RNA-capping assembly line. Nat Struct Mol Biol 14: 449-451. 

31. Egloflf MP, Benarroch D, Selisko B, RometteJL, Canard B (2002) An RNA cap 
(nucleoside-2'-0-)-methyltransferase in the flavivirus RNA polymerase NS5: 
crystal structure and functional characterization. EMBO J 21: 2757—2768. 

32. Egloflf MP, Decroly E, Malet H, Selisko B, Benarroch D, et al. (2007) Structural 
and functional analysis of methylation and 5'-RNA sequence requirements of 
short capped RNAs by the methyltransferase domain of dengue virus NS5. J Mol 
Biol 372: 723-736. 

33. Ray D, Shah A, Tilgner M, Guo Y, Zhao Y, et al. (2006) West Nile virus 5'-cap 
structure is formed by sequential guanine N-7 and ribose 2'-0 methylations by 
nonstructural protein 5. J Virol 80: 8362-8370. 

34. Bollati M, Alvarez K, Assenberg R, Baronti C, Canard B, et al. (2010) Structure 
and functionality in flavivirus NS-proteins: perspectives for drug design. 
Antiviral Res 87: 125-148. 

35. Snijder EJ, Bredenbeek PJ, Dobbe JC, Thiel V, Ziebuhr J, et al. (2003) Unique 
and conserved features of genome and proteome of SARS-coronavirus, an early 
split-off from the coronavirus group 2 lineage. J Mol Biol 331: 991—1004. 

36. Balzarini J, De Clercq E, Serafinowski P, Dorland E, Harrap KR (1992) 
Synthesis and antiviral activity of some new S-adenosyl-L-homocysteine 
derivatives. J Med Chem 35: 4576-4583. 

37. Dong H, Zhang B, Shi PY (2008) Flavivirus methyltransferase: a novel antiviral 
target. Antiviral Res 80: 1-10. 

38. Pugh CS, Borchardt RT, Stone HO (1977) Inhibition of Newcastle disease 
virion messenger RNA (guanine-7-)-methyltransferase by analogues of S- 
adenosylhomocysteine. Biochemistry 16: 3928—3932. 

39. Pugh CS, Borchardt RT, Stone HO (1978) Sinefungin, a potent inhibitor of 
virion mRNA(guanine-7-)-methyltransferase, mRNA(nucleoside-2'-(-methyl¬ 
transferase, and viral multiplication. J Biol Chem 253: 4075—4077. 

40. Debarnot C, Imbert I, Ferron F, Gluais L, Varlet I, et al. (2011) Crystallization 
and diffraction analysis of the SARS coronavirus nspl0/nspl6 complex. Acta 
Crystallogr Sect F Struct Biol Cryst Commun 67: 404-8. 

41. Holm L, Park J (2000) DaliLite workbench for protein structure comparison. 
Bioinformatics 16: 566-567. 

42. Dudev T, Lim C (2003) Principles governing Mg, Ca, and Zn binding and 
selectivity in proteins. Chem Rev 103: 773-788. 

43. Marcotrigiano J, Gingras AC, Sonenberg N, Burley SK (1997) Cocrystal 
structure of the messenger RNA 5' cap-binding protein (eIF4E) bound to 7- 
methyl-GDP. Cell 89: 951-961. 

44. Mazza C, Segref A, Mattaj IW, Cusack S (2002) Large-scale induced fit 
recognition of an m(7)GpppG cap analogue by the human nuclear cap-binding 
complex. EMBO J 21: 5548-5557. 


45. Guilligay D, Tarendeau F, Resa-Infante P, Coloma R, Crepin T, et al. (2008) 
The structural basis for cap binding by influenza virus polymerase subunit PB2. 
Nat Struct Mol Biol 15: 500-506. 

46. Quiocho FA, Hu G, Gershon PD (2000) Structural basis of mRNA cap 
recognition by proteins. Curr Opin Struct Biol 10: 78—86. 

47. Hodel AE, Gershon PD, Quiocho FA (1998) Structural basis for sequence- 
nonspecific recognition of 5'-capped mRNA by a cap-modifying enzyme. Mol 
Cell 1: 443-447. 

48. Li C, Xia Y, Gao X, Gershon PD (2004) Mechanism of RNA 2'-0-methylation: 
evidence that the catalytic lysine acts to steer rather than deprotonate the target 
nucleophile. Biochemistry 43: 5680-5687. 

49. Sigel RK, Sigel H (2010) A stability concept for metal ion coordination to single- 
stranded nucleic acids and affinities of individual sites. Acc Chem Res 43: 
974-984. 

50. Dey S, Pal A, Ghakrabarti P, Janin J (2010) The subunit interfaces of weakly 
associated homodimeric proteins. J Mol Biol 398: 146—160. 

51. Yoneyama M, Kikuchi M, Natsukawa T, Shinobu N, Imaizumi T, et al. (2004) 
The RNA helicase RIG-I has an essential function in double-stranded RNA- 
induced innate antiviral responses. Nat Immunol 5: 730—737. 

52. Malet H, Egloff MP, Selisko B, Butcher RE, Wright PJ, et al. (2007) Crystal 
structure of the RNA polymerase domain of the West Nile virus non-structural 
protein 5. J Biol Chem 282: 10678-10689. 

53. Dong H, Liu L, Zou G, Zhao Y, Li Z, et al. (2010) Structural and functional 
analyses of a conserved hydrophobic pocket of flavivirus methyltransferase. J Biol 
Chem 285: 32586-32595. 

54. McCoy AJ, Grosse-Kunstleve RW, Adams PD, Winn MD, Storoni LC, et al. 
(2007) Phaser crystallographic software. Journal of applied crystallography 40: 
658-674. 

55. Gowtan K (2010) Recent developments in classical density modification. Acta 
Crystallogr D Biol Crystallogr 66: 470-478. 

56. Cowtan K (2006) The Buccaneer software for automated model building. 1. 
Tracing protein chains. Acta Crystallogr D Biol Crystallogr 62: 1002—1011. 

57. Cohen SX, Morris PJ, Fernandez FJ, Ben Jelloul M, Kakaris M, et al. (2004) 
Towards complete validated models in the next generation of ARP/wARP. Acta 
Crystallogr D Biol Crystallogr 60: 2222-2229. 

58. Emsley P, Cowtan K (2004) Coot: model-building tools for molecular graphics. 
Acta crystallographica Section D, Biological crystallography 60: 2126—2132. 

59. Collaborative Computational Project N (1994) The CCP4 suite: programs for 
protein crystallography. Acta Crystallogr D Biol Crystallogr 50: 760—763. 

60. Saff EB, Kuijlaars ABJ (1997) Distributing Many Points on a Sphere. Math Intell 
19:5-11. 

61. Peyrane F, Selisko B, Decroly E, Vasseur JJ, Benarroch D, et al. (2007) High- 
yield production of short GpppA- and 7MeGpppA-capped RNAs and HPLC- 
monitoring of methyltransfer reactions at the guanine-N7 and adenosine-2'O 
positions. Nucleic Acids Res 35: e26. 



PLoS Pathogens | www.plospathogens.org 


14 


May 2011 


| Volume 7 | Issue 5 | el 002059 



Copyright of PLoS Pathogens is the property of Public Library of Science and its content may not be copied or 
emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. 
However, users may print, download, or email articles for individual use. 



